{
 "cells": [
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "# p21: Eigenvalues of Mathieu operator"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 5,
   "metadata": {
    "collapsed": false
   },
   "outputs": [],
   "source": [
    "%matplotlib inline\n",
    "%config InlineBackend.figure_format='svg'\n",
    "from numpy import pi,arange,sin,cos,zeros,diag,sort,real\n",
    "from scipy.linalg import toeplitz\n",
    "from numpy.linalg import eig\n",
    "from itertools import cycle\n",
    "from matplotlib.pyplot import figure,plot,xlabel,ylabel"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 6,
   "metadata": {
    "collapsed": false
   },
   "outputs": [
    {
     "data": {
      "image/svg+xml": [
       "<?xml version=\"1.0\" encoding=\"utf-8\" standalone=\"no\"?>\n",
       "<!DOCTYPE svg PUBLIC \"-//W3C//DTD SVG 1.1//EN\"\n",
       "  \"http://www.w3.org/Graphics/SVG/1.1/DTD/svg11.dtd\">\n",
       "<!-- Created with matplotlib (http://matplotlib.org/) -->\n",
       "<svg height=\"605pt\" version=\"1.1\" viewBox=\"0 0 339 605\" width=\"339pt\" xmlns=\"http://www.w3.org/2000/svg\" xmlns:xlink=\"http://www.w3.org/1999/xlink\">\n",
       " <defs>\n",
       "  <style type=\"text/css\">\n",
       "*{stroke-linecap:butt;stroke-linejoin:round;stroke-miterlimit:100000;}\n",
       "  </style>\n",
       " </defs>\n",
       " <g id=\"figure_1\">\n",
       "  <g id=\"patch_1\">\n",
       "   <path d=\"M -0 605.595312 \n",
       "L 339.545312 605.595312 \n",
       "L 339.545312 0 \n",
       "L -0 0 \n",
       "z\n",
       "\" style=\"fill:none;\"/>\n",
       "  </g>\n",
       "  <g id=\"axes_1\">\n",
       "   <g id=\"patch_2\">\n",
       "    <path d=\"M 46.982813 570.039062 \n",
       "L 325.982813 570.039062 \n",
       "L 325.982813 12.039062 \n",
       "L 46.982813 12.039062 \n",
       "z\n",
       "\" style=\"fill:#ffffff;\"/>\n",
       "   </g>\n",
       "   <g id=\"line2d_1\">\n",
       "    <path clip-path=\"url(#pb324bc428c)\" d=\"M 46.982813 330.896205 \n",
       "L 50.470312 331.054943 \n",
       "L 53.957813 331.523153 \n",
       "L 57.445312 332.27881 \n",
       "L 60.932813 333.290765 \n",
       "L 64.420312 334.52431 \n",
       "L 67.907813 335.945569 \n",
       "L 71.395313 337.524078 \n",
       "L 74.882812 339.233803 \n",
       "L 78.370312 341.053171 \n",
       "L 81.857812 342.964605 \n",
       "L 85.345313 344.953891 \n",
       "L 88.832813 347.00956 \n",
       "L 92.320312 349.122331 \n",
       "L 95.807813 351.284655 \n",
       "L 99.295312 353.490358 \n",
       "L 102.782813 355.734346 \n",
       "L 106.270313 358.01239 \n",
       "L 109.757812 360.320943 \n",
       "L 113.245313 362.657013 \n",
       "L 116.732812 365.018055 \n",
       "L 120.220313 367.401886 \n",
       "L 123.707813 369.80662 \n",
       "L 127.195312 372.230621 \n",
       "L 130.682813 374.672454 \n",
       "L 134.170312 377.130858 \n",
       "L 137.657813 379.604717 \n",
       "L 141.145313 382.093037 \n",
       "L 144.632812 384.594929 \n",
       "L 148.120313 387.109595 \n",
       "L 151.607812 389.636314 \n",
       "L 155.095313 392.174432 \n",
       "L 158.582813 394.723356 \n",
       "L 162.070312 397.282542 \n",
       "L 165.557813 399.851493 \n",
       "L 169.045312 402.429753 \n",
       "L 172.532813 405.0169 \n",
       "L 176.020312 407.612546 \n",
       "L 179.507812 410.216329 \n",
       "L 182.995313 412.827914 \n",
       "L 186.482812 415.44699 \n",
       "L 189.970313 418.073265 \n",
       "L 193.457812 420.706467 \n",
       "L 196.945312 423.346342 \n",
       "L 200.432813 425.992649 \n",
       "L 203.920312 428.645164 \n",
       "L 207.407813 431.303676 \n",
       "L 210.895312 433.967983 \n",
       "L 214.382813 436.637898 \n",
       "L 217.870313 439.313242 \n",
       "L 221.357812 441.993846 \n",
       "L 224.845313 444.679549 \n",
       "L 228.332812 447.3702 \n",
       "L 231.820313 450.065655 \n",
       "L 235.307813 452.765774 \n",
       "L 238.795312 455.470428 \n",
       "L 242.282813 458.179491 \n",
       "L 245.770312 460.892845 \n",
       "L 249.257813 463.610374 \n",
       "L 252.745313 466.331971 \n",
       "L 256.232813 469.057532 \n",
       "L 259.720313 471.786956 \n",
       "L 263.207812 474.520147 \n",
       "L 266.695313 477.257015 \n",
       "L 270.182813 479.997471 \n",
       "L 273.670313 482.741431 \n",
       "L 277.157813 485.488813 \n",
       "L 280.645312 488.239539 \n",
       "L 284.132813 490.993535 \n",
       "L 287.620313 493.750727 \n",
       "L 291.107813 496.511046 \n",
       "L 294.595313 499.274426 \n",
       "L 298.082812 502.0408 \n",
       "L 301.570313 504.810106 \n",
       "L 305.057813 507.582285 \n",
       "\" style=\"fill:none;stroke:#0000ff;stroke-linecap:square;\"/>\n",
       "   </g>\n",
       "   <g id=\"line2d_2\">\n",
       "    <path clip-path=\"url(#pb324bc428c)\" d=\"M 46.982813 322.924777 \n",
       "L 50.470312 324.557932 \n",
       "L 53.957813 326.26496 \n",
       "L 57.445312 328.040273 \n",
       "L 60.932813 329.878606 \n",
       "L 64.420312 331.775046 \n",
       "L 67.907813 333.725041 \n",
       "L 71.395313 335.724399 \n",
       "L 74.882812 337.76928 \n",
       "L 78.370312 339.856174 \n",
       "L 81.857812 341.981884 \n",
       "L 85.345313 344.1435 \n",
       "L 88.832813 346.33838 \n",
       "L 92.320312 348.564121 \n",
       "L 95.807813 350.818542 \n",
       "L 99.295312 353.099661 \n",
       "L 102.782813 355.405675 \n",
       "L 106.270313 357.734944 \n",
       "L 109.757812 360.085973 \n",
       "L 113.245313 362.4574 \n",
       "L 116.732812 364.847978 \n",
       "L 120.220313 367.256568 \n",
       "L 123.707813 369.682127 \n",
       "L 127.195312 372.123696 \n",
       "L 130.682813 374.580397 \n",
       "L 134.170312 377.051419 \n",
       "L 137.657813 379.536015 \n",
       "L 141.145313 382.033496 \n",
       "L 144.632812 384.543223 \n",
       "L 148.120313 387.064606 \n",
       "L 151.607812 389.597096 \n",
       "L 155.095313 392.140184 \n",
       "L 158.582813 394.693397 \n",
       "L 162.070312 397.256291 \n",
       "L 165.557813 399.828456 \n",
       "L 169.045312 402.409504 \n",
       "L 172.532813 404.999076 \n",
       "L 176.020312 407.596834 \n",
       "L 179.507812 410.202459 \n",
       "L 182.995313 412.815655 \n",
       "L 186.482812 415.43614 \n",
       "L 189.970313 418.06365 \n",
       "L 193.457812 420.697937 \n",
       "L 196.945312 423.338765 \n",
       "L 200.432813 425.985911 \n",
       "L 203.920312 428.639165 \n",
       "L 207.407813 431.298329 \n",
       "L 210.895312 433.963213 \n",
       "L 214.382813 436.633638 \n",
       "L 217.870313 439.309434 \n",
       "L 221.357812 441.990438 \n",
       "L 224.845313 444.676497 \n",
       "L 228.332812 447.367464 \n",
       "L 231.820313 450.0632 \n",
       "L 235.307813 452.763569 \n",
       "L 238.795312 455.468446 \n",
       "L 242.282813 458.177709 \n",
       "L 245.770312 460.89124 \n",
       "L 249.257813 463.608928 \n",
       "L 252.745313 466.330668 \n",
       "L 256.232813 469.056355 \n",
       "L 259.720313 471.785893 \n",
       "L 263.207812 474.519187 \n",
       "L 266.695313 477.256147 \n",
       "L 270.182813 479.996685 \n",
       "L 273.670313 482.740719 \n",
       "L 277.157813 485.488168 \n",
       "L 280.645312 488.238954 \n",
       "L 284.132813 490.993004 \n",
       "L 287.620313 493.750245 \n",
       "L 291.107813 496.510609 \n",
       "L 294.595313 499.274028 \n",
       "L 298.082812 502.040438 \n",
       "L 301.570313 504.809777 \n",
       "L 305.057813 507.581986 \n",
       "\" style=\"fill:none;stroke:#008000;stroke-dasharray:6.000000,6.000000;stroke-dashoffset:0.0;\"/>\n",
       "   </g>\n",
       "   <g id=\"line2d_3\">\n",
       "    <path clip-path=\"url(#pb324bc428c)\" d=\"M 46.982813 322.924777 \n",
       "L 50.470312 321.371352 \n",
       "L 53.957813 319.903711 \n",
       "L 57.445312 318.527993 \n",
       "L 60.932813 317.250287 \n",
       "L 64.420312 316.076458 \n",
       "L 67.907813 315.011962 \n",
       "L 71.395313 314.061632 \n",
       "L 74.882812 313.22948 \n",
       "L 78.370312 312.518511 \n",
       "L 81.857812 311.930583 \n",
       "L 85.345313 311.466318 \n",
       "L 88.832813 311.125083 \n",
       "L 92.320312 310.90504 \n",
       "L 95.807813 310.803256 \n",
       "L 99.295312 310.815865 \n",
       "L 102.782813 310.938256 \n",
       "L 106.270313 311.165275 \n",
       "L 109.757812 311.491418 \n",
       "L 113.245313 311.911014 \n",
       "L 116.732812 312.418369 \n",
       "L 120.220313 313.007893 \n",
       "L 123.707813 313.674188 \n",
       "L 127.195312 314.412113 \n",
       "L 130.682813 315.216824 \n",
       "L 134.170312 316.083796 \n",
       "L 137.657813 317.008832 \n",
       "L 141.145313 317.988056 \n",
       "L 144.632812 319.017905 \n",
       "L 148.120313 320.095108 \n",
       "L 151.607812 321.216674 \n",
       "L 155.095313 322.379864 \n",
       "L 158.582813 323.582178 \n",
       "L 162.070312 324.821326 \n",
       "L 165.557813 326.095218 \n",
       "L 169.045312 327.401937 \n",
       "L 172.532813 328.739729 \n",
       "L 176.020312 330.106982 \n",
       "L 179.507812 331.502219 \n",
       "L 182.995313 332.924076 \n",
       "L 186.482812 334.371299 \n",
       "L 189.970313 335.842728 \n",
       "L 193.457812 337.337291 \n",
       "L 196.945312 338.853995 \n",
       "L 200.432813 340.391916 \n",
       "L 203.920312 341.950198 \n",
       "L 207.407813 343.52804 \n",
       "L 210.895312 345.124697 \n",
       "L 214.382813 346.739471 \n",
       "L 217.870313 348.371709 \n",
       "L 221.357812 350.020798 \n",
       "L 224.845313 351.686162 \n",
       "L 228.332812 353.36726 \n",
       "L 231.820313 355.063581 \n",
       "L 235.307813 356.774644 \n",
       "L 238.795312 358.499992 \n",
       "L 242.282813 360.239195 \n",
       "L 245.770312 361.991845 \n",
       "L 249.257813 363.757553 \n",
       "L 252.745313 365.535953 \n",
       "L 256.232813 367.326693 \n",
       "L 259.720313 369.12944 \n",
       "L 263.207812 370.943877 \n",
       "L 266.695313 372.7697 \n",
       "L 270.182813 374.60662 \n",
       "L 273.670313 376.454359 \n",
       "L 277.157813 378.312653 \n",
       "L 280.645312 380.181248 \n",
       "L 284.132813 382.059901 \n",
       "L 287.620313 383.948378 \n",
       "L 291.107813 385.846455 \n",
       "L 294.595313 387.753917 \n",
       "L 298.082812 389.670557 \n",
       "L 301.570313 391.596175 \n",
       "L 305.057813 393.53058 \n",
       "\" style=\"fill:none;stroke:#ff0000;stroke-linecap:square;\"/>\n",
       "   </g>\n",
       "   <g id=\"line2d_4\">\n",
       "    <path clip-path=\"url(#pb324bc428c)\" d=\"M 46.982813 299.010491 \n",
       "L 50.470312 299.037058 \n",
       "L 53.957813 299.116703 \n",
       "L 57.445312 299.249262 \n",
       "L 60.932813 299.43446 \n",
       "L 64.420312 299.671922 \n",
       "L 67.907813 299.961169 \n",
       "L 71.395313 300.301628 \n",
       "L 74.882812 300.692637 \n",
       "L 78.370312 301.133455 \n",
       "L 81.857812 301.623265 \n",
       "L 85.345313 302.161183 \n",
       "L 88.832813 302.746271 \n",
       "L 92.320312 303.377538 \n",
       "L 95.807813 304.053954 \n",
       "L 99.295312 304.774456 \n",
       "L 102.782813 305.537954 \n",
       "L 106.270313 306.343342 \n",
       "L 109.757812 307.189501 \n",
       "L 113.245313 308.075308 \n",
       "L 116.732812 308.999639 \n",
       "L 120.220313 309.961379 \n",
       "L 123.707813 310.959418 \n",
       "L 127.195312 311.992666 \n",
       "L 130.682813 313.060047 \n",
       "L 134.170312 314.160506 \n",
       "L 137.657813 315.293012 \n",
       "L 141.145313 316.456558 \n",
       "L 144.632812 317.650163 \n",
       "L 148.120313 318.872874 \n",
       "L 151.607812 320.123767 \n",
       "L 155.095313 321.401946 \n",
       "L 158.582813 322.706544 \n",
       "L 162.070312 324.036725 \n",
       "L 165.557813 325.39168 \n",
       "L 169.045312 326.770629 \n",
       "L 172.532813 328.172822 \n",
       "L 176.020312 329.597536 \n",
       "L 179.507812 331.044075 \n",
       "L 182.995313 332.511769 \n",
       "L 186.482812 333.999975 \n",
       "L 189.970313 335.508074 \n",
       "L 193.457812 337.035473 \n",
       "L 196.945312 338.581599 \n",
       "L 200.432813 340.145905 \n",
       "L 203.920312 341.727863 \n",
       "L 207.407813 343.326968 \n",
       "L 210.895312 344.942733 \n",
       "L 214.382813 346.574691 \n",
       "L 217.870313 348.222393 \n",
       "L 221.357812 349.88541 \n",
       "L 224.845313 351.563325 \n",
       "L 228.332812 353.255742 \n",
       "L 231.820313 354.962278 \n",
       "L 235.307813 356.682564 \n",
       "L 238.795312 358.416248 \n",
       "L 242.282813 360.162987 \n",
       "L 245.770312 361.922455 \n",
       "L 249.257813 363.694337 \n",
       "L 252.745313 365.478329 \n",
       "L 256.232813 367.274138 \n",
       "L 259.720313 369.081483 \n",
       "L 263.207812 370.900092 \n",
       "L 266.695313 372.729704 \n",
       "L 270.182813 374.570066 \n",
       "L 273.670313 376.420934 \n",
       "L 277.157813 378.282075 \n",
       "L 280.645312 380.15326 \n",
       "L 284.132813 382.034272 \n",
       "L 287.620313 383.924898 \n",
       "L 291.107813 385.824934 \n",
       "L 294.595313 387.734182 \n",
       "L 298.082812 389.652452 \n",
       "L 301.570313 391.579558 \n",
       "L 305.057813 393.515321 \n",
       "\" style=\"fill:none;stroke:#00bfbf;stroke-dasharray:6.000000,6.000000;stroke-dashoffset:0.0;\"/>\n",
       "   </g>\n",
       "   <g id=\"line2d_5\">\n",
       "    <path clip-path=\"url(#pb324bc428c)\" d=\"M 46.982813 299.010491 \n",
       "L 50.470312 298.878332 \n",
       "L 53.957813 298.489933 \n",
       "L 57.445312 297.867552 \n",
       "L 60.932813 297.042723 \n",
       "L 64.420312 296.050692 \n",
       "L 67.907813 294.926019 \n",
       "L 71.395313 293.699997 \n",
       "L 74.882812 292.399629 \n",
       "L 78.370312 291.047587 \n",
       "L 81.857812 289.662675 \n",
       "L 85.345313 288.260439 \n",
       "L 88.832813 286.853786 \n",
       "L 92.320312 285.45352 \n",
       "L 95.807813 284.068772 \n",
       "L 99.295312 282.70735 \n",
       "L 102.782813 281.376001 \n",
       "L 106.270313 280.080605 \n",
       "L 109.757812 278.82632 \n",
       "L 113.245313 277.617692 \n",
       "L 116.732812 276.458723 \n",
       "L 120.220313 275.352926 \n",
       "L 123.707813 274.303356 \n",
       "L 127.195312 273.312637 \n",
       "L 130.682813 272.382973 \n",
       "L 134.170312 271.516159 \n",
       "L 137.657813 270.713595 \n",
       "L 141.145313 269.976292 \n",
       "L 144.632812 269.304887 \n",
       "L 148.120313 268.699659 \n",
       "L 151.607812 268.160549 \n",
       "L 155.095313 267.687182 \n",
       "L 158.582813 267.278899 \n",
       "L 162.070312 266.934781 \n",
       "L 165.557813 266.653682 \n",
       "L 169.045312 266.434267 \n",
       "L 172.532813 266.275037 \n",
       "L 176.020312 266.174367 \n",
       "L 179.507812 266.130533 \n",
       "L 182.995313 266.141743 \n",
       "L 186.482812 266.206159 \n",
       "L 189.970313 266.321921 \n",
       "L 193.457812 266.48717 \n",
       "L 196.945312 266.700061 \n",
       "L 200.432813 266.958776 \n",
       "L 203.920312 267.261541 \n",
       "L 207.407813 267.60663 \n",
       "L 210.895312 267.992371 \n",
       "L 214.382813 268.417156 \n",
       "L 217.870313 268.87944 \n",
       "L 221.357812 269.377743 \n",
       "L 224.845313 269.910653 \n",
       "L 228.332812 270.476824 \n",
       "L 231.820313 271.074976 \n",
       "L 235.307813 271.703893 \n",
       "L 238.795312 272.362419 \n",
       "L 242.282813 273.049461 \n",
       "L 245.770312 273.763981 \n",
       "L 249.257813 274.504997 \n",
       "L 252.745313 275.271579 \n",
       "L 256.232813 276.062846 \n",
       "L 259.720313 276.877964 \n",
       "L 263.207812 277.716142 \n",
       "L 266.695313 278.576633 \n",
       "L 270.182813 279.458727 \n",
       "L 273.670313 280.361751 \n",
       "L 277.157813 281.285069 \n",
       "L 280.645312 282.228073 \n",
       "L 284.132813 283.19019 \n",
       "L 287.620313 284.170872 \n",
       "L 291.107813 285.169601 \n",
       "L 294.595313 286.185882 \n",
       "L 298.082812 287.219243 \n",
       "L 301.570313 288.269236 \n",
       "L 305.057813 289.335433 \n",
       "\" style=\"fill:none;stroke:#bf00bf;stroke-linecap:square;\"/>\n",
       "   </g>\n",
       "   <g id=\"line2d_6\">\n",
       "    <path clip-path=\"url(#pb324bc428c)\" d=\"M 46.982813 259.153348 \n",
       "L 50.470312 259.134407 \n",
       "L 53.957813 259.081454 \n",
       "L 57.445312 259.000079 \n",
       "L 60.932813 258.895554 \n",
       "L 64.420312 258.772798 \n",
       "L 67.907813 258.636363 \n",
       "L 71.395313 258.490438 \n",
       "L 74.882812 258.338853 \n",
       "L 78.370312 258.185099 \n",
       "L 81.857812 258.032344 \n",
       "L 85.345313 257.883458 \n",
       "L 88.832813 257.741029 \n",
       "L 92.320312 257.607391 \n",
       "L 95.807813 257.484641 \n",
       "L 99.295312 257.374661 \n",
       "L 102.782813 257.279134 \n",
       "L 106.270313 257.199565 \n",
       "L 109.757812 257.137293 \n",
       "L 113.245313 257.093504 \n",
       "L 116.732812 257.069249 \n",
       "L 120.220313 257.065452 \n",
       "L 123.707813 257.082922 \n",
       "L 127.195312 257.122361 \n",
       "L 130.682813 257.184375 \n",
       "L 134.170312 257.269479 \n",
       "L 137.657813 257.378107 \n",
       "L 141.145313 257.510619 \n",
       "L 144.632812 257.667305 \n",
       "L 148.120313 257.848389 \n",
       "L 151.607812 258.054042 \n",
       "L 155.095313 258.284376 \n",
       "L 158.582813 258.539457 \n",
       "L 162.070312 258.819306 \n",
       "L 165.557813 259.123903 \n",
       "L 169.045312 259.45319 \n",
       "L 172.532813 259.807075 \n",
       "L 176.020312 260.185435 \n",
       "L 179.507812 260.588121 \n",
       "L 182.995313 261.014957 \n",
       "L 186.482812 261.465745 \n",
       "L 189.970313 261.940268 \n",
       "L 193.457812 262.43829 \n",
       "L 196.945312 262.959559 \n",
       "L 200.432813 263.503812 \n",
       "L 203.920312 264.07077 \n",
       "L 207.407813 264.660147 \n",
       "L 210.895312 265.271646 \n",
       "L 214.382813 265.904966 \n",
       "L 217.870313 266.559797 \n",
       "L 221.357812 267.235826 \n",
       "L 224.845313 267.932735 \n",
       "L 228.332812 268.650206 \n",
       "L 231.820313 269.387917 \n",
       "L 235.307813 270.145547 \n",
       "L 238.795312 270.922773 \n",
       "L 242.282813 271.719276 \n",
       "L 245.770312 272.534735 \n",
       "L 249.257813 273.368832 \n",
       "L 252.745313 274.221254 \n",
       "L 256.232813 275.091687 \n",
       "L 259.720313 275.979823 \n",
       "L 263.207812 276.885356 \n",
       "L 266.695313 277.807985 \n",
       "L 270.182813 278.747412 \n",
       "L 273.670313 279.703346 \n",
       "L 277.157813 280.675496 \n",
       "L 280.645312 281.66358 \n",
       "L 284.132813 282.66732 \n",
       "L 287.620313 283.68644 \n",
       "L 291.107813 284.720673 \n",
       "L 294.595313 285.769755 \n",
       "L 298.082812 286.833427 \n",
       "L 301.570313 287.911435 \n",
       "L 305.057813 289.003532 \n",
       "\" style=\"fill:none;stroke:#bfbf00;stroke-dasharray:6.000000,6.000000;stroke-dashoffset:0.0;\"/>\n",
       "   </g>\n",
       "   <g id=\"line2d_7\">\n",
       "    <path clip-path=\"url(#pb324bc428c)\" d=\"M 46.982813 259.153348 \n",
       "L 50.470312 259.132416 \n",
       "L 53.957813 259.065561 \n",
       "L 57.445312 258.946653 \n",
       "L 60.932813 258.769624 \n",
       "L 64.420312 258.528637 \n",
       "L 67.907813 258.218282 \n",
       "L 71.395313 257.833783 \n",
       "L 74.882812 257.371207 \n",
       "L 78.370312 256.827642 \n",
       "L 81.857812 256.201349 \n",
       "L 85.345313 255.491846 \n",
       "L 88.832813 254.69993 \n",
       "L 92.320312 253.827631 \n",
       "L 95.807813 252.878101 \n",
       "L 99.295312 251.855455 \n",
       "L 102.782813 250.764585 \n",
       "L 106.270313 249.610957 \n",
       "L 109.757812 248.400419 \n",
       "L 113.245313 247.139023 \n",
       "L 116.732812 245.832875 \n",
       "L 120.220313 244.488014 \n",
       "L 123.707813 243.110322 \n",
       "L 127.195312 241.705459 \n",
       "L 130.682813 240.27882 \n",
       "L 134.170312 238.835516 \n",
       "L 137.657813 237.380362 \n",
       "L 141.145313 235.917883 \n",
       "L 144.632812 234.452319 \n",
       "L 148.120313 232.987639 \n",
       "L 151.607812 231.52756 \n",
       "L 155.095313 230.07556 \n",
       "L 158.582813 228.634898 \n",
       "L 162.070312 227.208624 \n",
       "L 165.557813 225.7996 \n",
       "L 169.045312 224.410504 \n",
       "L 172.532813 223.043849 \n",
       "L 176.020312 221.701984 \n",
       "L 179.507812 220.387107 \n",
       "L 182.995313 219.101264 \n",
       "L 186.482812 217.846359 \n",
       "L 189.970313 216.624151 \n",
       "L 193.457812 215.436261 \n",
       "L 196.945312 214.284166 \n",
       "L 200.432813 213.169208 \n",
       "L 203.920312 212.092587 \n",
       "L 207.407813 211.055365 \n",
       "L 210.895312 210.058465 \n",
       "L 214.382813 209.102671 \n",
       "L 217.870313 208.188634 \n",
       "L 221.357812 207.316867 \n",
       "L 224.845313 206.487755 \n",
       "L 228.332812 205.701554 \n",
       "L 231.820313 204.958396 \n",
       "L 235.307813 204.258297 \n",
       "L 238.795312 203.601162 \n",
       "L 242.282813 202.98679 \n",
       "L 245.770312 202.414882 \n",
       "L 249.257813 201.885053 \n",
       "L 252.745313 201.396834 \n",
       "L 256.232813 200.949684 \n",
       "L 259.720313 200.542997 \n",
       "L 263.207812 200.176112 \n",
       "L 266.695313 199.84832 \n",
       "L 270.182813 199.55887 \n",
       "L 273.670313 199.306981 \n",
       "L 277.157813 199.091844 \n",
       "L 280.645312 198.912633 \n",
       "L 284.132813 198.768508 \n",
       "L 287.620313 198.658622 \n",
       "L 291.107813 198.582126 \n",
       "L 294.595313 198.538171 \n",
       "L 298.082812 198.525915 \n",
       "L 301.570313 198.544526 \n",
       "L 305.057813 198.593184 \n",
       "\" style=\"fill:none;stroke:#000000;stroke-linecap:square;\"/>\n",
       "   </g>\n",
       "   <g id=\"line2d_8\">\n",
       "    <path clip-path=\"url(#pb324bc428c)\" d=\"M 46.982813 203.353348 \n",
       "L 50.470312 203.342724 \n",
       "L 53.957813 203.310909 \n",
       "L 57.445312 203.258069 \n",
       "L 60.932813 203.184481 \n",
       "L 64.420312 203.09053 \n",
       "L 67.907813 202.976698 \n",
       "L 71.395313 202.843567 \n",
       "L 74.882812 202.691809 \n",
       "L 78.370312 202.522174 \n",
       "L 81.857812 202.335493 \n",
       "L 85.345313 202.132658 \n",
       "L 88.832813 201.914623 \n",
       "L 92.320312 201.682389 \n",
       "L 95.807813 201.436999 \n",
       "L 99.295312 201.17953 \n",
       "L 102.782813 200.911082 \n",
       "L 106.270313 200.632774 \n",
       "L 109.757812 200.345733 \n",
       "L 113.245313 200.051092 \n",
       "L 116.732812 199.749981 \n",
       "L 120.220313 199.443521 \n",
       "L 123.707813 199.132823 \n",
       "L 127.195312 198.818977 \n",
       "L 130.682813 198.503057 \n",
       "L 134.170312 198.186109 \n",
       "L 137.657813 197.869156 \n",
       "L 141.145313 197.55319 \n",
       "L 144.632812 197.239172 \n",
       "L 148.120313 196.928032 \n",
       "L 151.607812 196.620667 \n",
       "L 155.095313 196.317937 \n",
       "L 158.582813 196.02067 \n",
       "L 162.070312 195.729659 \n",
       "L 165.557813 195.44566 \n",
       "L 169.045312 195.169395 \n",
       "L 172.532813 194.901551 \n",
       "L 176.020312 194.642779 \n",
       "L 179.507812 194.393698 \n",
       "L 182.995313 194.154892 \n",
       "L 186.482812 193.926913 \n",
       "L 189.970313 193.710279 \n",
       "L 193.457812 193.505478 \n",
       "L 196.945312 193.312966 \n",
       "L 200.432813 193.133171 \n",
       "L 203.920312 192.966489 \n",
       "L 207.407813 192.813291 \n",
       "L 210.895312 192.673917 \n",
       "L 214.382813 192.548684 \n",
       "L 217.870313 192.43788 \n",
       "L 221.357812 192.341771 \n",
       "L 224.845313 192.260597 \n",
       "L 228.332812 192.194577 \n",
       "L 231.820313 192.143907 \n",
       "L 235.307813 192.10876 \n",
       "L 238.795312 192.089293 \n",
       "L 242.282813 192.085638 \n",
       "L 245.770312 192.097914 \n",
       "L 249.257813 192.126217 \n",
       "L 252.745313 192.170629 \n",
       "L 256.232813 192.231215 \n",
       "L 259.720313 192.308025 \n",
       "L 263.207812 192.401095 \n",
       "L 266.695313 192.510444 \n",
       "L 270.182813 192.636081 \n",
       "L 273.670313 192.778002 \n",
       "L 277.157813 192.93619 \n",
       "L 280.645312 193.110617 \n",
       "L 284.132813 193.301247 \n",
       "L 287.620313 193.508032 \n",
       "L 291.107813 193.730915 \n",
       "L 294.595313 193.96983 \n",
       "L 298.082812 194.224705 \n",
       "L 301.570313 194.49546 \n",
       "L 305.057813 194.782006 \n",
       "\" style=\"fill:none;stroke:#0000ff;stroke-dasharray:6.000000,6.000000;stroke-dashoffset:0.0;\"/>\n",
       "   </g>\n",
       "   <g id=\"line2d_9\">\n",
       "    <path clip-path=\"url(#pb324bc428c)\" d=\"M 46.982813 203.353348 \n",
       "L 50.470312 203.342713 \n",
       "L 53.957813 203.310732 \n",
       "L 57.445312 203.257174 \n",
       "L 60.932813 203.181659 \n",
       "L 64.420312 203.083656 \n",
       "L 67.907813 202.962487 \n",
       "L 71.395313 202.817333 \n",
       "L 74.882812 202.647237 \n",
       "L 78.370312 202.451113 \n",
       "L 81.857812 202.227752 \n",
       "L 85.345313 201.975838 \n",
       "L 88.832813 201.693957 \n",
       "L 92.320312 201.380612 \n",
       "L 95.807813 201.034246 \n",
       "L 99.295312 200.65326 \n",
       "L 102.782813 200.236036 \n",
       "L 106.270313 199.780967 \n",
       "L 109.757812 199.286482 \n",
       "L 113.245313 198.751081 \n",
       "L 116.732812 198.173363 \n",
       "L 120.220313 197.552062 \n",
       "L 123.707813 196.886073 \n",
       "L 127.195312 196.174489 \n",
       "L 130.682813 195.416621 \n",
       "L 134.170312 194.612026 \n",
       "L 137.657813 193.760518 \n",
       "L 141.145313 192.862189 \n",
       "L 144.632812 191.917407 \n",
       "L 148.120313 190.926817 \n",
       "L 151.607812 189.891336 \n",
       "L 155.095313 188.812136 \n",
       "L 158.582813 187.690632 \n",
       "L 162.070312 186.528454 \n",
       "L 165.557813 185.327422 \n",
       "L 169.045312 184.089525 \n",
       "L 172.532813 182.816885 \n",
       "L 176.020312 181.511735 \n",
       "L 179.507812 180.176391 \n",
       "L 182.995313 178.813224 \n",
       "L 186.482812 177.424642 \n",
       "L 189.970313 176.013066 \n",
       "L 193.457812 174.580914 \n",
       "L 196.945312 173.130584 \n",
       "L 200.432813 171.664444 \n",
       "L 203.920312 170.184821 \n",
       "L 207.407813 168.69399 \n",
       "L 210.895312 167.194175 \n",
       "L 214.382813 165.687537 \n",
       "L 217.870313 164.176175 \n",
       "L 221.357812 162.662125 \n",
       "L 224.845313 161.147358 \n",
       "L 228.332812 159.633777 \n",
       "L 231.820313 158.123224 \n",
       "L 235.307813 156.617474 \n",
       "L 238.795312 155.118241 \n",
       "L 242.282813 153.627176 \n",
       "L 245.770312 152.145871 \n",
       "L 249.257813 150.675856 \n",
       "L 252.745313 149.218604 \n",
       "L 256.232813 147.775532 \n",
       "L 259.720313 146.347997 \n",
       "L 263.207812 144.937302 \n",
       "L 266.695313 143.544693 \n",
       "L 270.182813 142.17136 \n",
       "L 273.670313 140.818438 \n",
       "L 277.157813 139.487006 \n",
       "L 280.645312 138.178084 \n",
       "L 284.132813 136.892638 \n",
       "L 287.620313 135.631576 \n",
       "L 291.107813 134.395748 \n",
       "L 294.595313 133.185943 \n",
       "L 298.082812 132.002896 \n",
       "L 301.570313 130.847277 \n",
       "L 305.057813 129.719701 \n",
       "\" style=\"fill:none;stroke:#008000;stroke-linecap:square;\"/>\n",
       "   </g>\n",
       "   <g id=\"line2d_10\">\n",
       "    <path clip-path=\"url(#pb324bc428c)\" d=\"M 46.982813 131.610491 \n",
       "L 50.470312 131.603848 \n",
       "L 53.957813 131.583917 \n",
       "L 57.445312 131.550695 \n",
       "L 60.932813 131.504177 \n",
       "L 64.420312 131.44436 \n",
       "L 67.907813 131.371246 \n",
       "L 71.395313 131.284842 \n",
       "L 74.882812 131.185163 \n",
       "L 78.370312 131.072231 \n",
       "L 81.857812 130.94608 \n",
       "L 85.345313 130.806758 \n",
       "L 88.832813 130.654325 \n",
       "L 92.320312 130.488857 \n",
       "L 95.807813 130.310445 \n",
       "L 99.295312 130.119198 \n",
       "L 102.782813 129.915244 \n",
       "L 106.270313 129.698728 \n",
       "L 109.757812 129.469818 \n",
       "L 113.245313 129.228699 \n",
       "L 116.732812 128.975576 \n",
       "L 120.220313 128.710677 \n",
       "L 123.707813 128.434249 \n",
       "L 127.195312 128.146558 \n",
       "L 130.682813 127.847892 \n",
       "L 134.170312 127.538557 \n",
       "L 137.657813 127.21888 \n",
       "L 141.145313 126.889204 \n",
       "L 144.632812 126.549889 \n",
       "L 148.120313 126.201315 \n",
       "L 151.607812 125.843875 \n",
       "L 155.095313 125.477976 \n",
       "L 158.582813 125.104041 \n",
       "L 162.070312 124.722504 \n",
       "L 165.557813 124.333809 \n",
       "L 169.045312 123.938412 \n",
       "L 172.532813 123.536778 \n",
       "L 176.020312 123.129377 \n",
       "L 179.507812 122.716688 \n",
       "L 182.995313 122.299193 \n",
       "L 186.482812 121.877381 \n",
       "L 189.970313 121.45174 \n",
       "L 193.457812 121.022763 \n",
       "L 196.945312 120.590941 \n",
       "L 200.432813 120.156766 \n",
       "L 203.920312 119.720729 \n",
       "L 207.407813 119.283317 \n",
       "L 210.895312 118.845015 \n",
       "L 214.382813 118.406304 \n",
       "L 217.870313 117.967658 \n",
       "L 221.357812 117.52955 \n",
       "L 224.845313 117.09244 \n",
       "L 228.332812 116.656788 \n",
       "L 231.820313 116.223041 \n",
       "L 235.307813 115.791641 \n",
       "L 238.795312 115.363021 \n",
       "L 242.282813 114.937604 \n",
       "L 245.770312 114.515804 \n",
       "L 249.257813 114.098026 \n",
       "L 252.745313 113.684666 \n",
       "L 256.232813 113.276106 \n",
       "L 259.720313 112.872722 \n",
       "L 263.207812 112.474877 \n",
       "L 266.695313 112.082922 \n",
       "L 270.182813 111.697201 \n",
       "L 273.670313 111.318042 \n",
       "L 277.157813 110.945766 \n",
       "L 280.645312 110.580682 \n",
       "L 284.132813 110.223086 \n",
       "L 287.620313 109.873265 \n",
       "L 291.107813 109.531494 \n",
       "L 294.595313 109.198036 \n",
       "L 298.082812 108.873146 \n",
       "L 301.570313 108.557064 \n",
       "L 305.057813 108.250023 \n",
       "\" style=\"fill:none;stroke:#ff0000;stroke-dasharray:6.000000,6.000000;stroke-dashoffset:0.0;\"/>\n",
       "   </g>\n",
       "   <g id=\"line2d_11\">\n",
       "    <path clip-path=\"url(#pb324bc428c)\" d=\"M 46.982813 131.610491 \n",
       "L 50.470312 131.603848 \n",
       "L 53.957813 131.583916 \n",
       "L 57.445312 131.550686 \n",
       "L 60.932813 131.504141 \n",
       "L 64.420312 131.444252 \n",
       "L 67.907813 131.370978 \n",
       "L 71.395313 131.284264 \n",
       "L 74.882812 131.184038 \n",
       "L 78.370312 131.070208 \n",
       "L 81.857812 130.942662 \n",
       "L 85.345313 130.801268 \n",
       "L 88.832813 130.645867 \n",
       "L 92.320312 130.476275 \n",
       "L 95.807813 130.29228 \n",
       "L 99.295312 130.093644 \n",
       "L 102.782813 129.880094 \n",
       "L 106.270313 129.651331 \n",
       "L 109.757812 129.407019 \n",
       "L 113.245313 129.146793 \n",
       "L 116.732812 128.870252 \n",
       "L 120.220313 128.576961 \n",
       "L 123.707813 128.266451 \n",
       "L 127.195312 127.938222 \n",
       "L 130.682813 127.591738 \n",
       "L 134.170312 127.226431 \n",
       "L 137.657813 126.841703 \n",
       "L 141.145313 126.436927 \n",
       "L 144.632812 126.011451 \n",
       "L 148.120313 125.564596 \n",
       "L 151.607812 125.095666 \n",
       "L 155.095313 124.603945 \n",
       "L 158.582813 124.088707 \n",
       "L 162.070312 123.54922 \n",
       "L 165.557813 122.984748 \n",
       "L 169.045312 122.394561 \n",
       "L 172.532813 121.777941 \n",
       "L 176.020312 121.134186 \n",
       "L 179.507812 120.462622 \n",
       "L 182.995313 119.762608 \n",
       "L 186.482812 119.033543 \n",
       "L 189.970313 118.274877 \n",
       "L 193.457812 117.486116 \n",
       "L 196.945312 116.66683 \n",
       "L 200.432813 115.816661 \n",
       "L 203.920312 114.93533 \n",
       "L 207.407813 114.022641 \n",
       "L 210.895312 113.078488 \n",
       "L 214.382813 112.102857 \n",
       "L 217.870313 111.095831 \n",
       "L 221.357812 110.057592 \n",
       "L 224.845313 108.988417 \n",
       "L 228.332812 107.888686 \n",
       "L 231.820313 106.758872 \n",
       "L 235.307813 105.599543 \n",
       "L 238.795312 104.411358 \n",
       "L 242.282813 103.19506 \n",
       "L 245.770312 101.951472 \n",
       "L 249.257813 100.681493 \n",
       "L 252.745313 99.386087 \n",
       "L 256.232813 98.066281 \n",
       "L 259.720313 96.723155 \n",
       "L 263.207812 95.357833 \n",
       "L 266.695313 93.971484 \n",
       "L 270.182813 92.565305 \n",
       "L 273.670313 91.140521 \n",
       "L 277.157813 89.698379 \n",
       "L 280.645312 88.240136 \n",
       "L 284.132813 86.767063 \n",
       "L 287.620313 85.280431 \n",
       "L 291.107813 83.781514 \n",
       "L 294.595313 82.27158 \n",
       "L 298.082812 80.75189 \n",
       "L 301.570313 79.223693 \n",
       "L 305.057813 77.688228 \n",
       "\" style=\"fill:none;stroke:#00bfbf;stroke-linecap:square;\"/>\n",
       "   </g>\n",
       "   <g id=\"patch_3\">\n",
       "    <path d=\"M 46.982813 12.039062 \n",
       "L 325.982813 12.039062 \n",
       "\" style=\"fill:none;stroke:#000000;stroke-linecap:square;stroke-linejoin:miter;\"/>\n",
       "   </g>\n",
       "   <g id=\"patch_4\">\n",
       "    <path d=\"M 325.982813 570.039062 \n",
       "L 325.982813 12.039062 \n",
       "\" style=\"fill:none;stroke:#000000;stroke-linecap:square;stroke-linejoin:miter;\"/>\n",
       "   </g>\n",
       "   <g id=\"patch_5\">\n",
       "    <path d=\"M 46.982813 570.039062 \n",
       "L 325.982813 570.039062 \n",
       "\" style=\"fill:none;stroke:#000000;stroke-linecap:square;stroke-linejoin:miter;\"/>\n",
       "   </g>\n",
       "   <g id=\"patch_6\">\n",
       "    <path d=\"M 46.982813 570.039062 \n",
       "L 46.982813 12.039062 \n",
       "\" style=\"fill:none;stroke:#000000;stroke-linecap:square;stroke-linejoin:miter;\"/>\n",
       "   </g>\n",
       "   <g id=\"matplotlib.axis_1\">\n",
       "    <g id=\"xtick_1\">\n",
       "     <g id=\"line2d_12\">\n",
       "      <defs>\n",
       "       <path d=\"M 0 0 \n",
       "L 0 -4 \n",
       "\" id=\"m9bcfe46652\" style=\"stroke:#000000;stroke-width:0.5;\"/>\n",
       "      </defs>\n",
       "      <g>\n",
       "       <use style=\"stroke:#000000;stroke-width:0.5;\" x=\"46.9828125\" xlink:href=\"#m9bcfe46652\" y=\"570.0390625\"/>\n",
       "      </g>\n",
       "     </g>\n",
       "     <g id=\"line2d_13\">\n",
       "      <defs>\n",
       "       <path d=\"M 0 0 \n",
       "L 0 4 \n",
       "\" id=\"m0d4d8ed4dd\" style=\"stroke:#000000;stroke-width:0.5;\"/>\n",
       "      </defs>\n",
       "      <g>\n",
       "       <use style=\"stroke:#000000;stroke-width:0.5;\" x=\"46.9828125\" xlink:href=\"#m0d4d8ed4dd\" y=\"12.0390625\"/>\n",
       "      </g>\n",
       "     </g>\n",
       "     <g id=\"text_1\">\n",
       "      <!-- 0 -->\n",
       "      <defs>\n",
       "       <path d=\"M 31.78125 66.40625 \n",
       "Q 24.171875 66.40625 20.328125 58.90625 \n",
       "Q 16.5 51.421875 16.5 36.375 \n",
       "Q 16.5 21.390625 20.328125 13.890625 \n",
       "Q 24.171875 6.390625 31.78125 6.390625 \n",
       "Q 39.453125 6.390625 43.28125 13.890625 \n",
       "Q 47.125 21.390625 47.125 36.375 \n",
       "Q 47.125 51.421875 43.28125 58.90625 \n",
       "Q 39.453125 66.40625 31.78125 66.40625 \n",
       "M 31.78125 74.21875 \n",
       "Q 44.046875 74.21875 50.515625 64.515625 \n",
       "Q 56.984375 54.828125 56.984375 36.375 \n",
       "Q 56.984375 17.96875 50.515625 8.265625 \n",
       "Q 44.046875 -1.421875 31.78125 -1.421875 \n",
       "Q 19.53125 -1.421875 13.0625 8.265625 \n",
       "Q 6.59375 17.96875 6.59375 36.375 \n",
       "Q 6.59375 54.828125 13.0625 64.515625 \n",
       "Q 19.53125 74.21875 31.78125 74.21875 \n",
       "\" id=\"BitstreamVeraSans-Roman-30\"/>\n",
       "      </defs>\n",
       "      <g transform=\"translate(43.8015625 581.6375)scale(0.1 -0.1)\">\n",
       "       <use xlink:href=\"#BitstreamVeraSans-Roman-30\"/>\n",
       "      </g>\n",
       "     </g>\n",
       "    </g>\n",
       "    <g id=\"xtick_2\">\n",
       "     <g id=\"line2d_14\">\n",
       "      <g>\n",
       "       <use style=\"stroke:#000000;stroke-width:0.5;\" x=\"81.8578125\" xlink:href=\"#m9bcfe46652\" y=\"570.0390625\"/>\n",
       "      </g>\n",
       "     </g>\n",
       "     <g id=\"line2d_15\">\n",
       "      <g>\n",
       "       <use style=\"stroke:#000000;stroke-width:0.5;\" x=\"81.8578125\" xlink:href=\"#m0d4d8ed4dd\" y=\"12.0390625\"/>\n",
       "      </g>\n",
       "     </g>\n",
       "     <g id=\"text_2\">\n",
       "      <!-- 2 -->\n",
       "      <defs>\n",
       "       <path d=\"M 19.1875 8.296875 \n",
       "L 53.609375 8.296875 \n",
       "L 53.609375 0 \n",
       "L 7.328125 0 \n",
       "L 7.328125 8.296875 \n",
       "Q 12.9375 14.109375 22.625 23.890625 \n",
       "Q 32.328125 33.6875 34.8125 36.53125 \n",
       "Q 39.546875 41.84375 41.421875 45.53125 \n",
       "Q 43.3125 49.21875 43.3125 52.78125 \n",
       "Q 43.3125 58.59375 39.234375 62.25 \n",
       "Q 35.15625 65.921875 28.609375 65.921875 \n",
       "Q 23.96875 65.921875 18.8125 64.3125 \n",
       "Q 13.671875 62.703125 7.8125 59.421875 \n",
       "L 7.8125 69.390625 \n",
       "Q 13.765625 71.78125 18.9375 73 \n",
       "Q 24.125 74.21875 28.421875 74.21875 \n",
       "Q 39.75 74.21875 46.484375 68.546875 \n",
       "Q 53.21875 62.890625 53.21875 53.421875 \n",
       "Q 53.21875 48.921875 51.53125 44.890625 \n",
       "Q 49.859375 40.875 45.40625 35.40625 \n",
       "Q 44.1875 33.984375 37.640625 27.21875 \n",
       "Q 31.109375 20.453125 19.1875 8.296875 \n",
       "\" id=\"BitstreamVeraSans-Roman-32\"/>\n",
       "      </defs>\n",
       "      <g transform=\"translate(78.6765625 581.6375)scale(0.1 -0.1)\">\n",
       "       <use xlink:href=\"#BitstreamVeraSans-Roman-32\"/>\n",
       "      </g>\n",
       "     </g>\n",
       "    </g>\n",
       "    <g id=\"xtick_3\">\n",
       "     <g id=\"line2d_16\">\n",
       "      <g>\n",
       "       <use style=\"stroke:#000000;stroke-width:0.5;\" x=\"116.7328125\" xlink:href=\"#m9bcfe46652\" y=\"570.0390625\"/>\n",
       "      </g>\n",
       "     </g>\n",
       "     <g id=\"line2d_17\">\n",
       "      <g>\n",
       "       <use style=\"stroke:#000000;stroke-width:0.5;\" x=\"116.7328125\" xlink:href=\"#m0d4d8ed4dd\" y=\"12.0390625\"/>\n",
       "      </g>\n",
       "     </g>\n",
       "     <g id=\"text_3\">\n",
       "      <!-- 4 -->\n",
       "      <defs>\n",
       "       <path d=\"M 37.796875 64.3125 \n",
       "L 12.890625 25.390625 \n",
       "L 37.796875 25.390625 \n",
       "z\n",
       "M 35.203125 72.90625 \n",
       "L 47.609375 72.90625 \n",
       "L 47.609375 25.390625 \n",
       "L 58.015625 25.390625 \n",
       "L 58.015625 17.1875 \n",
       "L 47.609375 17.1875 \n",
       "L 47.609375 0 \n",
       "L 37.796875 0 \n",
       "L 37.796875 17.1875 \n",
       "L 4.890625 17.1875 \n",
       "L 4.890625 26.703125 \n",
       "z\n",
       "\" id=\"BitstreamVeraSans-Roman-34\"/>\n",
       "      </defs>\n",
       "      <g transform=\"translate(113.5515625 581.6375)scale(0.1 -0.1)\">\n",
       "       <use xlink:href=\"#BitstreamVeraSans-Roman-34\"/>\n",
       "      </g>\n",
       "     </g>\n",
       "    </g>\n",
       "    <g id=\"xtick_4\">\n",
       "     <g id=\"line2d_18\">\n",
       "      <g>\n",
       "       <use style=\"stroke:#000000;stroke-width:0.5;\" x=\"151.6078125\" xlink:href=\"#m9bcfe46652\" y=\"570.0390625\"/>\n",
       "      </g>\n",
       "     </g>\n",
       "     <g id=\"line2d_19\">\n",
       "      <g>\n",
       "       <use style=\"stroke:#000000;stroke-width:0.5;\" x=\"151.6078125\" xlink:href=\"#m0d4d8ed4dd\" y=\"12.0390625\"/>\n",
       "      </g>\n",
       "     </g>\n",
       "     <g id=\"text_4\">\n",
       "      <!-- 6 -->\n",
       "      <defs>\n",
       "       <path d=\"M 33.015625 40.375 \n",
       "Q 26.375 40.375 22.484375 35.828125 \n",
       "Q 18.609375 31.296875 18.609375 23.390625 \n",
       "Q 18.609375 15.53125 22.484375 10.953125 \n",
       "Q 26.375 6.390625 33.015625 6.390625 \n",
       "Q 39.65625 6.390625 43.53125 10.953125 \n",
       "Q 47.40625 15.53125 47.40625 23.390625 \n",
       "Q 47.40625 31.296875 43.53125 35.828125 \n",
       "Q 39.65625 40.375 33.015625 40.375 \n",
       "M 52.59375 71.296875 \n",
       "L 52.59375 62.3125 \n",
       "Q 48.875 64.0625 45.09375 64.984375 \n",
       "Q 41.3125 65.921875 37.59375 65.921875 \n",
       "Q 27.828125 65.921875 22.671875 59.328125 \n",
       "Q 17.53125 52.734375 16.796875 39.40625 \n",
       "Q 19.671875 43.65625 24.015625 45.921875 \n",
       "Q 28.375 48.1875 33.59375 48.1875 \n",
       "Q 44.578125 48.1875 50.953125 41.515625 \n",
       "Q 57.328125 34.859375 57.328125 23.390625 \n",
       "Q 57.328125 12.15625 50.6875 5.359375 \n",
       "Q 44.046875 -1.421875 33.015625 -1.421875 \n",
       "Q 20.359375 -1.421875 13.671875 8.265625 \n",
       "Q 6.984375 17.96875 6.984375 36.375 \n",
       "Q 6.984375 53.65625 15.1875 63.9375 \n",
       "Q 23.390625 74.21875 37.203125 74.21875 \n",
       "Q 40.921875 74.21875 44.703125 73.484375 \n",
       "Q 48.484375 72.75 52.59375 71.296875 \n",
       "\" id=\"BitstreamVeraSans-Roman-36\"/>\n",
       "      </defs>\n",
       "      <g transform=\"translate(148.4265625 581.6375)scale(0.1 -0.1)\">\n",
       "       <use xlink:href=\"#BitstreamVeraSans-Roman-36\"/>\n",
       "      </g>\n",
       "     </g>\n",
       "    </g>\n",
       "    <g id=\"xtick_5\">\n",
       "     <g id=\"line2d_20\">\n",
       "      <g>\n",
       "       <use style=\"stroke:#000000;stroke-width:0.5;\" x=\"186.4828125\" xlink:href=\"#m9bcfe46652\" y=\"570.0390625\"/>\n",
       "      </g>\n",
       "     </g>\n",
       "     <g id=\"line2d_21\">\n",
       "      <g>\n",
       "       <use style=\"stroke:#000000;stroke-width:0.5;\" x=\"186.4828125\" xlink:href=\"#m0d4d8ed4dd\" y=\"12.0390625\"/>\n",
       "      </g>\n",
       "     </g>\n",
       "     <g id=\"text_5\">\n",
       "      <!-- 8 -->\n",
       "      <defs>\n",
       "       <path d=\"M 31.78125 34.625 \n",
       "Q 24.75 34.625 20.71875 30.859375 \n",
       "Q 16.703125 27.09375 16.703125 20.515625 \n",
       "Q 16.703125 13.921875 20.71875 10.15625 \n",
       "Q 24.75 6.390625 31.78125 6.390625 \n",
       "Q 38.8125 6.390625 42.859375 10.171875 \n",
       "Q 46.921875 13.96875 46.921875 20.515625 \n",
       "Q 46.921875 27.09375 42.890625 30.859375 \n",
       "Q 38.875 34.625 31.78125 34.625 \n",
       "M 21.921875 38.8125 \n",
       "Q 15.578125 40.375 12.03125 44.71875 \n",
       "Q 8.5 49.078125 8.5 55.328125 \n",
       "Q 8.5 64.0625 14.71875 69.140625 \n",
       "Q 20.953125 74.21875 31.78125 74.21875 \n",
       "Q 42.671875 74.21875 48.875 69.140625 \n",
       "Q 55.078125 64.0625 55.078125 55.328125 \n",
       "Q 55.078125 49.078125 51.53125 44.71875 \n",
       "Q 48 40.375 41.703125 38.8125 \n",
       "Q 48.828125 37.15625 52.796875 32.3125 \n",
       "Q 56.78125 27.484375 56.78125 20.515625 \n",
       "Q 56.78125 9.90625 50.3125 4.234375 \n",
       "Q 43.84375 -1.421875 31.78125 -1.421875 \n",
       "Q 19.734375 -1.421875 13.25 4.234375 \n",
       "Q 6.78125 9.90625 6.78125 20.515625 \n",
       "Q 6.78125 27.484375 10.78125 32.3125 \n",
       "Q 14.796875 37.15625 21.921875 38.8125 \n",
       "M 18.3125 54.390625 \n",
       "Q 18.3125 48.734375 21.84375 45.5625 \n",
       "Q 25.390625 42.390625 31.78125 42.390625 \n",
       "Q 38.140625 42.390625 41.71875 45.5625 \n",
       "Q 45.3125 48.734375 45.3125 54.390625 \n",
       "Q 45.3125 60.0625 41.71875 63.234375 \n",
       "Q 38.140625 66.40625 31.78125 66.40625 \n",
       "Q 25.390625 66.40625 21.84375 63.234375 \n",
       "Q 18.3125 60.0625 18.3125 54.390625 \n",
       "\" id=\"BitstreamVeraSans-Roman-38\"/>\n",
       "      </defs>\n",
       "      <g transform=\"translate(183.3015625 581.6375)scale(0.1 -0.1)\">\n",
       "       <use xlink:href=\"#BitstreamVeraSans-Roman-38\"/>\n",
       "      </g>\n",
       "     </g>\n",
       "    </g>\n",
       "    <g id=\"xtick_6\">\n",
       "     <g id=\"line2d_22\">\n",
       "      <g>\n",
       "       <use style=\"stroke:#000000;stroke-width:0.5;\" x=\"221.3578125\" xlink:href=\"#m9bcfe46652\" y=\"570.0390625\"/>\n",
       "      </g>\n",
       "     </g>\n",
       "     <g id=\"line2d_23\">\n",
       "      <g>\n",
       "       <use style=\"stroke:#000000;stroke-width:0.5;\" x=\"221.3578125\" xlink:href=\"#m0d4d8ed4dd\" y=\"12.0390625\"/>\n",
       "      </g>\n",
       "     </g>\n",
       "     <g id=\"text_6\">\n",
       "      <!-- 10 -->\n",
       "      <defs>\n",
       "       <path d=\"M 12.40625 8.296875 \n",
       "L 28.515625 8.296875 \n",
       "L 28.515625 63.921875 \n",
       "L 10.984375 60.40625 \n",
       "L 10.984375 69.390625 \n",
       "L 28.421875 72.90625 \n",
       "L 38.28125 72.90625 \n",
       "L 38.28125 8.296875 \n",
       "L 54.390625 8.296875 \n",
       "L 54.390625 0 \n",
       "L 12.40625 0 \n",
       "z\n",
       "\" id=\"BitstreamVeraSans-Roman-31\"/>\n",
       "      </defs>\n",
       "      <g transform=\"translate(214.9953125 581.6375)scale(0.1 -0.1)\">\n",
       "       <use xlink:href=\"#BitstreamVeraSans-Roman-31\"/>\n",
       "       <use x=\"63.623046875\" xlink:href=\"#BitstreamVeraSans-Roman-30\"/>\n",
       "      </g>\n",
       "     </g>\n",
       "    </g>\n",
       "    <g id=\"xtick_7\">\n",
       "     <g id=\"line2d_24\">\n",
       "      <g>\n",
       "       <use style=\"stroke:#000000;stroke-width:0.5;\" x=\"256.2328125\" xlink:href=\"#m9bcfe46652\" y=\"570.0390625\"/>\n",
       "      </g>\n",
       "     </g>\n",
       "     <g id=\"line2d_25\">\n",
       "      <g>\n",
       "       <use style=\"stroke:#000000;stroke-width:0.5;\" x=\"256.2328125\" xlink:href=\"#m0d4d8ed4dd\" y=\"12.0390625\"/>\n",
       "      </g>\n",
       "     </g>\n",
       "     <g id=\"text_7\">\n",
       "      <!-- 12 -->\n",
       "      <g transform=\"translate(249.8703125 581.6375)scale(0.1 -0.1)\">\n",
       "       <use xlink:href=\"#BitstreamVeraSans-Roman-31\"/>\n",
       "       <use x=\"63.623046875\" xlink:href=\"#BitstreamVeraSans-Roman-32\"/>\n",
       "      </g>\n",
       "     </g>\n",
       "    </g>\n",
       "    <g id=\"xtick_8\">\n",
       "     <g id=\"line2d_26\">\n",
       "      <g>\n",
       "       <use style=\"stroke:#000000;stroke-width:0.5;\" x=\"291.1078125\" xlink:href=\"#m9bcfe46652\" y=\"570.0390625\"/>\n",
       "      </g>\n",
       "     </g>\n",
       "     <g id=\"line2d_27\">\n",
       "      <g>\n",
       "       <use style=\"stroke:#000000;stroke-width:0.5;\" x=\"291.1078125\" xlink:href=\"#m0d4d8ed4dd\" y=\"12.0390625\"/>\n",
       "      </g>\n",
       "     </g>\n",
       "     <g id=\"text_8\">\n",
       "      <!-- 14 -->\n",
       "      <g transform=\"translate(284.7453125 581.6375)scale(0.1 -0.1)\">\n",
       "       <use xlink:href=\"#BitstreamVeraSans-Roman-31\"/>\n",
       "       <use x=\"63.623046875\" xlink:href=\"#BitstreamVeraSans-Roman-34\"/>\n",
       "      </g>\n",
       "     </g>\n",
       "    </g>\n",
       "    <g id=\"xtick_9\">\n",
       "     <g id=\"line2d_28\">\n",
       "      <g>\n",
       "       <use style=\"stroke:#000000;stroke-width:0.5;\" x=\"325.9828125\" xlink:href=\"#m9bcfe46652\" y=\"570.0390625\"/>\n",
       "      </g>\n",
       "     </g>\n",
       "     <g id=\"line2d_29\">\n",
       "      <g>\n",
       "       <use style=\"stroke:#000000;stroke-width:0.5;\" x=\"325.9828125\" xlink:href=\"#m0d4d8ed4dd\" y=\"12.0390625\"/>\n",
       "      </g>\n",
       "     </g>\n",
       "     <g id=\"text_9\">\n",
       "      <!-- 16 -->\n",
       "      <g transform=\"translate(319.6203125 581.6375)scale(0.1 -0.1)\">\n",
       "       <use xlink:href=\"#BitstreamVeraSans-Roman-31\"/>\n",
       "       <use x=\"63.623046875\" xlink:href=\"#BitstreamVeraSans-Roman-36\"/>\n",
       "      </g>\n",
       "     </g>\n",
       "    </g>\n",
       "    <g id=\"text_10\">\n",
       "     <!-- q -->\n",
       "     <defs>\n",
       "      <path d=\"M 14.796875 27.296875 \n",
       "Q 14.796875 17.390625 18.875 11.75 \n",
       "Q 22.953125 6.109375 30.078125 6.109375 \n",
       "Q 37.203125 6.109375 41.296875 11.75 \n",
       "Q 45.40625 17.390625 45.40625 27.296875 \n",
       "Q 45.40625 37.203125 41.296875 42.84375 \n",
       "Q 37.203125 48.484375 30.078125 48.484375 \n",
       "Q 22.953125 48.484375 18.875 42.84375 \n",
       "Q 14.796875 37.203125 14.796875 27.296875 \n",
       "M 45.40625 8.203125 \n",
       "Q 42.578125 3.328125 38.25 0.953125 \n",
       "Q 33.9375 -1.421875 27.875 -1.421875 \n",
       "Q 17.96875 -1.421875 11.734375 6.484375 \n",
       "Q 5.515625 14.40625 5.515625 27.296875 \n",
       "Q 5.515625 40.1875 11.734375 48.09375 \n",
       "Q 17.96875 56 27.875 56 \n",
       "Q 33.9375 56 38.25 53.625 \n",
       "Q 42.578125 51.265625 45.40625 46.390625 \n",
       "L 45.40625 54.6875 \n",
       "L 54.390625 54.6875 \n",
       "L 54.390625 -20.796875 \n",
       "L 45.40625 -20.796875 \n",
       "z\n",
       "\" id=\"BitstreamVeraSans-Roman-71\"/>\n",
       "     </defs>\n",
       "     <g transform=\"translate(183.30859375 596.315625)scale(0.1 -0.1)\">\n",
       "      <use xlink:href=\"#BitstreamVeraSans-Roman-71\"/>\n",
       "     </g>\n",
       "    </g>\n",
       "   </g>\n",
       "   <g id=\"matplotlib.axis_2\">\n",
       "    <g id=\"ytick_1\">\n",
       "     <g id=\"line2d_30\">\n",
       "      <defs>\n",
       "       <path d=\"M 0 0 \n",
       "L 4 0 \n",
       "\" id=\"m2d31cba656\" style=\"stroke:#000000;stroke-width:0.5;\"/>\n",
       "      </defs>\n",
       "      <g>\n",
       "       <use style=\"stroke:#000000;stroke-width:0.5;\" x=\"46.9828125\" xlink:href=\"#m2d31cba656\" y=\"570.0390625\"/>\n",
       "      </g>\n",
       "     </g>\n",
       "     <g id=\"line2d_31\">\n",
       "      <defs>\n",
       "       <path d=\"M 0 0 \n",
       "L -4 0 \n",
       "\" id=\"m90faafee00\" style=\"stroke:#000000;stroke-width:0.5;\"/>\n",
       "      </defs>\n",
       "      <g>\n",
       "       <use style=\"stroke:#000000;stroke-width:0.5;\" x=\"325.9828125\" xlink:href=\"#m90faafee00\" y=\"570.0390625\"/>\n",
       "      </g>\n",
       "     </g>\n",
       "     <g id=\"text_11\">\n",
       "      <!-- −30 -->\n",
       "      <defs>\n",
       "       <path d=\"M 10.59375 35.5 \n",
       "L 73.1875 35.5 \n",
       "L 73.1875 27.203125 \n",
       "L 10.59375 27.203125 \n",
       "z\n",
       "\" id=\"BitstreamVeraSans-Roman-2212\"/>\n",
       "       <path d=\"M 40.578125 39.3125 \n",
       "Q 47.65625 37.796875 51.625 33 \n",
       "Q 55.609375 28.21875 55.609375 21.1875 \n",
       "Q 55.609375 10.40625 48.1875 4.484375 \n",
       "Q 40.765625 -1.421875 27.09375 -1.421875 \n",
       "Q 22.515625 -1.421875 17.65625 -0.515625 \n",
       "Q 12.796875 0.390625 7.625 2.203125 \n",
       "L 7.625 11.71875 \n",
       "Q 11.71875 9.328125 16.59375 8.109375 \n",
       "Q 21.484375 6.890625 26.8125 6.890625 \n",
       "Q 36.078125 6.890625 40.9375 10.546875 \n",
       "Q 45.796875 14.203125 45.796875 21.1875 \n",
       "Q 45.796875 27.640625 41.28125 31.265625 \n",
       "Q 36.765625 34.90625 28.71875 34.90625 \n",
       "L 20.21875 34.90625 \n",
       "L 20.21875 43.015625 \n",
       "L 29.109375 43.015625 \n",
       "Q 36.375 43.015625 40.234375 45.921875 \n",
       "Q 44.09375 48.828125 44.09375 54.296875 \n",
       "Q 44.09375 59.90625 40.109375 62.90625 \n",
       "Q 36.140625 65.921875 28.71875 65.921875 \n",
       "Q 24.65625 65.921875 20.015625 65.03125 \n",
       "Q 15.375 64.15625 9.8125 62.3125 \n",
       "L 9.8125 71.09375 \n",
       "Q 15.4375 72.65625 20.34375 73.4375 \n",
       "Q 25.25 74.21875 29.59375 74.21875 \n",
       "Q 40.828125 74.21875 47.359375 69.109375 \n",
       "Q 53.90625 64.015625 53.90625 55.328125 \n",
       "Q 53.90625 49.265625 50.4375 45.09375 \n",
       "Q 46.96875 40.921875 40.578125 39.3125 \n",
       "\" id=\"BitstreamVeraSans-Roman-33\"/>\n",
       "      </defs>\n",
       "      <g transform=\"translate(21.878125 572.7984375)scale(0.1 -0.1)\">\n",
       "       <use xlink:href=\"#BitstreamVeraSans-Roman-2212\"/>\n",
       "       <use x=\"83.7890625\" xlink:href=\"#BitstreamVeraSans-Roman-33\"/>\n",
       "       <use x=\"147.412109375\" xlink:href=\"#BitstreamVeraSans-Roman-30\"/>\n",
       "      </g>\n",
       "     </g>\n",
       "    </g>\n",
       "    <g id=\"ytick_2\">\n",
       "     <g id=\"line2d_32\">\n",
       "      <g>\n",
       "       <use style=\"stroke:#000000;stroke-width:0.5;\" x=\"46.9828125\" xlink:href=\"#m2d31cba656\" y=\"490.324776786\"/>\n",
       "      </g>\n",
       "     </g>\n",
       "     <g id=\"line2d_33\">\n",
       "      <g>\n",
       "       <use style=\"stroke:#000000;stroke-width:0.5;\" x=\"325.9828125\" xlink:href=\"#m90faafee00\" y=\"490.324776786\"/>\n",
       "      </g>\n",
       "     </g>\n",
       "     <g id=\"text_12\">\n",
       "      <!-- −20 -->\n",
       "      <g transform=\"translate(21.878125 493.084151786)scale(0.1 -0.1)\">\n",
       "       <use xlink:href=\"#BitstreamVeraSans-Roman-2212\"/>\n",
       "       <use x=\"83.7890625\" xlink:href=\"#BitstreamVeraSans-Roman-32\"/>\n",
       "       <use x=\"147.412109375\" xlink:href=\"#BitstreamVeraSans-Roman-30\"/>\n",
       "      </g>\n",
       "     </g>\n",
       "    </g>\n",
       "    <g id=\"ytick_3\">\n",
       "     <g id=\"line2d_34\">\n",
       "      <g>\n",
       "       <use style=\"stroke:#000000;stroke-width:0.5;\" x=\"46.9828125\" xlink:href=\"#m2d31cba656\" y=\"410.610491071\"/>\n",
       "      </g>\n",
       "     </g>\n",
       "     <g id=\"line2d_35\">\n",
       "      <g>\n",
       "       <use style=\"stroke:#000000;stroke-width:0.5;\" x=\"325.9828125\" xlink:href=\"#m90faafee00\" y=\"410.610491071\"/>\n",
       "      </g>\n",
       "     </g>\n",
       "     <g id=\"text_13\">\n",
       "      <!-- −10 -->\n",
       "      <g transform=\"translate(21.878125 413.369866071)scale(0.1 -0.1)\">\n",
       "       <use xlink:href=\"#BitstreamVeraSans-Roman-2212\"/>\n",
       "       <use x=\"83.7890625\" xlink:href=\"#BitstreamVeraSans-Roman-31\"/>\n",
       "       <use x=\"147.412109375\" xlink:href=\"#BitstreamVeraSans-Roman-30\"/>\n",
       "      </g>\n",
       "     </g>\n",
       "    </g>\n",
       "    <g id=\"ytick_4\">\n",
       "     <g id=\"line2d_36\">\n",
       "      <g>\n",
       "       <use style=\"stroke:#000000;stroke-width:0.5;\" x=\"46.9828125\" xlink:href=\"#m2d31cba656\" y=\"330.896205357\"/>\n",
       "      </g>\n",
       "     </g>\n",
       "     <g id=\"line2d_37\">\n",
       "      <g>\n",
       "       <use style=\"stroke:#000000;stroke-width:0.5;\" x=\"325.9828125\" xlink:href=\"#m90faafee00\" y=\"330.896205357\"/>\n",
       "      </g>\n",
       "     </g>\n",
       "     <g id=\"text_14\">\n",
       "      <!-- 0 -->\n",
       "      <g transform=\"translate(36.6203125 333.655580357)scale(0.1 -0.1)\">\n",
       "       <use xlink:href=\"#BitstreamVeraSans-Roman-30\"/>\n",
       "      </g>\n",
       "     </g>\n",
       "    </g>\n",
       "    <g id=\"ytick_5\">\n",
       "     <g id=\"line2d_38\">\n",
       "      <g>\n",
       "       <use style=\"stroke:#000000;stroke-width:0.5;\" x=\"46.9828125\" xlink:href=\"#m2d31cba656\" y=\"251.181919643\"/>\n",
       "      </g>\n",
       "     </g>\n",
       "     <g id=\"line2d_39\">\n",
       "      <g>\n",
       "       <use style=\"stroke:#000000;stroke-width:0.5;\" x=\"325.9828125\" xlink:href=\"#m90faafee00\" y=\"251.181919643\"/>\n",
       "      </g>\n",
       "     </g>\n",
       "     <g id=\"text_15\">\n",
       "      <!-- 10 -->\n",
       "      <g transform=\"translate(30.2578125 253.941294643)scale(0.1 -0.1)\">\n",
       "       <use xlink:href=\"#BitstreamVeraSans-Roman-31\"/>\n",
       "       <use x=\"63.623046875\" xlink:href=\"#BitstreamVeraSans-Roman-30\"/>\n",
       "      </g>\n",
       "     </g>\n",
       "    </g>\n",
       "    <g id=\"ytick_6\">\n",
       "     <g id=\"line2d_40\">\n",
       "      <g>\n",
       "       <use style=\"stroke:#000000;stroke-width:0.5;\" x=\"46.9828125\" xlink:href=\"#m2d31cba656\" y=\"171.467633929\"/>\n",
       "      </g>\n",
       "     </g>\n",
       "     <g id=\"line2d_41\">\n",
       "      <g>\n",
       "       <use style=\"stroke:#000000;stroke-width:0.5;\" x=\"325.9828125\" xlink:href=\"#m90faafee00\" y=\"171.467633929\"/>\n",
       "      </g>\n",
       "     </g>\n",
       "     <g id=\"text_16\">\n",
       "      <!-- 20 -->\n",
       "      <g transform=\"translate(30.2578125 174.227008929)scale(0.1 -0.1)\">\n",
       "       <use xlink:href=\"#BitstreamVeraSans-Roman-32\"/>\n",
       "       <use x=\"63.623046875\" xlink:href=\"#BitstreamVeraSans-Roman-30\"/>\n",
       "      </g>\n",
       "     </g>\n",
       "    </g>\n",
       "    <g id=\"ytick_7\">\n",
       "     <g id=\"line2d_42\">\n",
       "      <g>\n",
       "       <use style=\"stroke:#000000;stroke-width:0.5;\" x=\"46.9828125\" xlink:href=\"#m2d31cba656\" y=\"91.7533482143\"/>\n",
       "      </g>\n",
       "     </g>\n",
       "     <g id=\"line2d_43\">\n",
       "      <g>\n",
       "       <use style=\"stroke:#000000;stroke-width:0.5;\" x=\"325.9828125\" xlink:href=\"#m90faafee00\" y=\"91.7533482143\"/>\n",
       "      </g>\n",
       "     </g>\n",
       "     <g id=\"text_17\">\n",
       "      <!-- 30 -->\n",
       "      <g transform=\"translate(30.2578125 94.5127232143)scale(0.1 -0.1)\">\n",
       "       <use xlink:href=\"#BitstreamVeraSans-Roman-33\"/>\n",
       "       <use x=\"63.623046875\" xlink:href=\"#BitstreamVeraSans-Roman-30\"/>\n",
       "      </g>\n",
       "     </g>\n",
       "    </g>\n",
       "    <g id=\"ytick_8\">\n",
       "     <g id=\"line2d_44\">\n",
       "      <g>\n",
       "       <use style=\"stroke:#000000;stroke-width:0.5;\" x=\"46.9828125\" xlink:href=\"#m2d31cba656\" y=\"12.0390625\"/>\n",
       "      </g>\n",
       "     </g>\n",
       "     <g id=\"line2d_45\">\n",
       "      <g>\n",
       "       <use style=\"stroke:#000000;stroke-width:0.5;\" x=\"325.9828125\" xlink:href=\"#m90faafee00\" y=\"12.0390625\"/>\n",
       "      </g>\n",
       "     </g>\n",
       "     <g id=\"text_18\">\n",
       "      <!-- 40 -->\n",
       "      <g transform=\"translate(30.2578125 14.7984375)scale(0.1 -0.1)\">\n",
       "       <use xlink:href=\"#BitstreamVeraSans-Roman-34\"/>\n",
       "       <use x=\"63.623046875\" xlink:href=\"#BitstreamVeraSans-Roman-30\"/>\n",
       "      </g>\n",
       "     </g>\n",
       "    </g>\n",
       "    <g id=\"text_19\">\n",
       "     <!-- $\\lambda$ -->\n",
       "     <defs>\n",
       "      <path d=\"M 5.078125 1.703125 \n",
       "Q 5.078125 3.171875 6.390625 4.5 \n",
       "L 34.515625 32.078125 \n",
       "L 24.609375 59.90625 \n",
       "Q 23.4375 63.09375 22.140625 64.9375 \n",
       "Q 20.84375 66.796875 18.40625 66.796875 \n",
       "Q 17.28125 66.796875 17.28125 68.109375 \n",
       "Q 17.390625 68.703125 17.78125 69.046875 \n",
       "Q 18.171875 69.390625 18.796875 69.390625 \n",
       "Q 29.296875 69.390625 31.5 63.1875 \n",
       "L 51.8125 6.5 \n",
       "Q 53.125 2.25 54.5 0.78125 \n",
       "Q 54.984375 0.046875 54.984375 -0.203125 \n",
       "Q 54.984375 -1.125 54.203125 -1.125 \n",
       "L 48.78125 -1.125 \n",
       "Q 46.296875 -0.390625 44.828125 3.21875 \n",
       "L 35.6875 28.8125 \n",
       "L 11.53125 0.390625 \n",
       "Q 9.71875 -1.3125 8.296875 -1.3125 \n",
       "Q 7.03125 -1.3125 6.046875 -0.4375 \n",
       "Q 5.078125 0.4375 5.078125 1.703125 \n",
       "\" id=\"Cmmi10-b8\"/>\n",
       "     </defs>\n",
       "     <g transform=\"translate(14.7984375 293.9890625)rotate(-90.0)scale(0.1 -0.1)\">\n",
       "      <use transform=\"translate(0.0 0.609375)\" xlink:href=\"#Cmmi10-b8\"/>\n",
       "     </g>\n",
       "    </g>\n",
       "   </g>\n",
       "  </g>\n",
       " </g>\n",
       " <defs>\n",
       "  <clipPath id=\"pb324bc428c\">\n",
       "   <rect height=\"558.0\" width=\"279.0\" x=\"46.9828125\" y=\"12.0390625\"/>\n",
       "  </clipPath>\n",
       " </defs>\n",
       "</svg>\n"
      ],
      "text/plain": [
       "<matplotlib.figure.Figure at 0x106aa8190>"
      ]
     },
     "metadata": {},
     "output_type": "display_data"
    }
   ],
   "source": [
    "N = 42; h = 2.0*pi/N; x = h*arange(1,N+1)\n",
    "col = zeros(N)\n",
    "col[0] = -pi**2/(3.0*h**2) - 1.0/6.0\n",
    "col[1:] = -0.5*(-1.0)**arange(1,N)/sin(0.5*h*arange(1,N))**2\n",
    "D2 = toeplitz(col)\n",
    "\n",
    "ne = 11 # number of eigenvalues to plot\n",
    "qq = arange(0.0, 15.0, 0.2)\n",
    "data= zeros((len(qq),ne))\n",
    "i = 0\n",
    "for q in qq:\n",
    "    evals,evecs = eig(-D2 + 2.0*q*diag(cos(2.0*x)))\n",
    "    e = real(sort(evals))\n",
    "    data[i,:] = e[0:ne]\n",
    "    i = i + 1\n",
    "    \n",
    "figure(figsize=(5,10))\n",
    "lines=cycle([\"-\",\"--\"])\n",
    "for i in range(ne):\n",
    "    plot(qq,data[:,i],next(lines))\n",
    "xlabel(\"q\")\n",
    "ylabel(\"$\\lambda$\");"
   ]
  }
 ],
 "metadata": {
  "anaconda-cloud": {},
  "kernelspec": {
   "display_name": "Python [default]",
   "language": "python",
   "name": "python2"
  },
  "language_info": {
   "codemirror_mode": {
    "name": "ipython",
    "version": 2
   },
   "file_extension": ".py",
   "mimetype": "text/x-python",
   "name": "python",
   "nbconvert_exporter": "python",
   "pygments_lexer": "ipython2",
   "version": "2.7.12"
  }
 },
 "nbformat": 4,
 "nbformat_minor": 0
}
